#include <iostream>
#include "interpolation_polynomial.h"
#include <vector>
using namespace std;

int main()
{
  vector<double> x = {0,6,10,13,17,20,28};
  vector<double> y = {6.67,17.3,42.7,37.3,30.1,29.3,28.7};
  Newton_formula p(x,y);
  vector<vector<double>> f = p.getf();
  for(int i = 0; i < x.size();i++)
    {
      cout<< f[i][i]<< " ";
      cout<< endl;
    }
      cout<< endl;
  return 0;
}
